Ballistic and diffuse transport through a ferromagnetic domain wall 
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We study transport through ballistic and diffuse ferromag- 
netic domain walls in a two-band Stoner model with a rotat- 
ing magnetization direction. For a ballistic domain wall, the 
change in the conductance due to the domain wall scatter- 
ing is obtained from an adiabatic approximation valid when 
the length of the domain wall is much longer than the Fermi 
wavelength. In diffuse systems, the change in the resistivity is 
calculated using a diagrammatic technique to the lowest order 
in the domain wall scattering and taking into account spin- 
dependent scattering lifetimes and screening of the domain 
wall potential. 
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I. INTRODUCTION 

In a ferroniagnet, domains with different directions of 
the magnetization are favored by the long-range magnetic 
dipolar interaction. The boundary between the domains, 
the domain walls (DW's), are a source of magnetoresis- 
tance that recently^ has attracted experimentalETEl and 
theoretical interest .BiS 

For ballistic systems, where the electron mean free 
path is longer than the system size, first-principles band- 
structure calculations have shown that the DW resis- 
tance is enhanced ckie to the nearly degenerate bands 
at the Fermi energy.S The rotating magnetization direc- 
tion causes an effective potential barrier for the electrons 
which increase the resistance. Recently large DW resis- 
tance in ballistic Ni nanocontacts has been measurecu in 
agreement with the predictions in Ref. 0. _. 

In the diffuse transport regime Cabrera and FalicovEHl 
interpreted transport through a single DW as a tunneling 
process and the corresponding MR was found to be ex- 
ponentially small. Bergerliil modelled the domain wall 
scattering as a force on the magnetic moment of the 
conduction electrons. Tatara and FukuyamatI calculated 
the DW conductivity for spin-independent scattering life- 
times. Levy and ZhangQ pointed out that spin-dependent 



impurity scattering can strongly enhance the DW resis- 
tivity. 

Beyond the .semiclassical transport theories, Tatara 
and FukuyamaD predict a negative DW resistivity as a 
result of the reduced weak localization correction due to 
the decoherence of the electrons by the scattering off the 
domain wall. However, quantum interference effects do 
not explain the experiments in Refs. bH, where the nega- 
tive DW resistance persists up to high temperature where 
the inelastic scattering length is shorter than the-^nean 
free path. It has been suggested by Ruediger et al.t3 that 
the experimentally observed negative DW resistance is an 
extrinsic effect caused by the interplay between orbital 
effects due to the internal magnetic fields and surface 
scattering. Recently, it was also demonstrated that the 
large negative domain wall resistance of Co filmsEI is due 
to MR resistivity anisotropyJla 

It is the purpose of the present paper to give a de- 
tailed account of the transport through a domain wall 
both in the ballistic and the diffuse regime in a two- 
band Stoner model. In the ballistic transport regime, 
the transport through the magnetic domain wall can be 
treated by an adiabatic approximation similar to the one 
used for transport through a quantum point contact. In 
the diffuse regime we will use the diagrammatic tech- 
nique introduced in Ref. @ and generalize it to the case 
of asymmetric impiirity scattering life-times (but without 
localization effecta3) and screening of the domain wall 
potential. Our results, although more general, reduce in 
the case of strong spin splitting to results that are very 
similar to those obtained in Ref. H using a Boltzmann 
equation. We explain why the results of the two meth- 
ods differ. Some results have been publishejdi-already in 
a brief report and a conference proceedings .Bo Here we 
give an in depth discussion of the results including the 
technical details of the derivations. 

The paper is organized in the following way. The two- 
band Stoner model for the ferromagnet with a rotating 
magnetic field and how it can be reduced to a more 
tractable form by a local gauge transformation is dis- 
cussed in section nil T he ballistic transport regime is 
discussed in sectionplll and the diffuse transport regime 
in section IV . We give our conclusions in section M. The 



appendices include the adiabatic approximation which 
can be used in the baUistic situation, recipes for the cal- 
culation of the frequency summation of the Feynman di- 
agrams, and the spin-spiral case which can be exactly 
diagonalized. 



this regime, the transport properties can be described by 
the Landauer conductance, 
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II. MODEL 

Throughout, we will use an effective 2-band model to 
describe the ferromagnet with the Hamiltonian 



H= f dr^lir) 



2m 



*o(r), (1) 



where hb is the Bohr magneton, a^ , <7y and cr^ are the 
Pauli spin matrices, H(r) is the effective magnetic field 
arising from the electronic exchange interaction, the mag- 
netic dipole interaction, and the external magnetic field 
and ^0 (r) is the 2-component spinor wave function. The 
direction of the effective magnetic field is represented by 
a rotation angle 9{z) which varies along the z-direction. 
The spin-orbit interaction and the Lorentz force due to 
the internal magnetization are disregarded, since exper- 
imentally the DW magnetoresistance can be separated 
from the anomalous magnetoresistance (AMR) and the 
orbital magnetp|resistance (OMR) □ We use a local gauge 
transformationQ '4'o(r) = [/(r)'J'(r), where 



where T^'^ is the transmission probability for an elec- 
tron in the transverse mode k|| to pass the domain wall 
from spin-state s' to spin-state s. Domain walls in tran- 
sition metals are much thicker than the Fermi wave- 
length (the length of the domain wall is A^, = 40nm 
(Fe), A^, — lOOnm (Ni), A^ = 15nm (Co) and the Fermi 
wavelength is roughly Xp ~ 0.2 nm). The transmis- 
sion probability can therefore be calculated with the aid 
of an adiabatic approximation on the eigenstates of the 
Hamiltonian after the gauge transformation (g) (see Ap- 
pendix H). The domain wall is then equivalent to an 
effective potential barrier for the electrons. The conduc- 
tance is determined by the minimum number of prop- 
agating modes, which is where the gradient of the ro- 
tating magnetic field has its maximum value, a{z) = 
dO{z)/dz -^ flniax (for details see Appendix]^. The con- 
ductance can then be found from the conductance of a 
spin-spiral with gradient flmax, which has the dispersion 

the conductance is 



n^l{2mW + <ax/4 ± (fc^aLx + P^Y'^ and 



U{y) = cos(6'/2)ct^ + sin(6i/2)cr^ 
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and introduce the Fourier transform of the change of 
the direction of the magnetic field a{z) = d6{z)/dz = 



^ exp(igz)aq. After the gauge transformation (g) the 

Hamiltonian (|l|) becomes H = U^HU = Hq + V. The 
unperturbed term is 



^o^Yl ('^k - A*) cLcks , 



(3) 



ks 



where s = + (s = — ) denotes spin-up (spin-down) states, 
e^ = ?i k'^ /2m~ sA and the spiBpSplitting is A = /iB|H|. 
The interaction with the DW ia3 
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where z is a unit vector in the z-direction and k, = kz. 



III. BALLISTIC TRANSPORT 

Transport through the domain wall is ballistic when 
the system size is smaller than the inean free path. In 



where Vg = dEs/ {Tidk^) is the group velocity. Carry- 
ing out the integration, we find the domain wall resis- 
tance_ i?^/i?o = R - Rq {l/R = G, 1/Rq = Go = 

{e^Akl)/{27rh)): 



Rw 

Rq 



P^^/kp -p^/{a- 



max i^/ 



max 
in ax 



<2p2 
> 2p2 



(7) 



where h^kf,/{2m) = Ep and A ee h^p^/{2m). The 
screening of the domain wall potential discussed below 
is not important for the calculation of the conductance, 
since by a calculation following the lines in section IV, we 



find that the shift in the chemical potential due to the 
rotating magnetization is 6^ ~ — (l/48)i?tt,(A/i?i?)^ -t- 
0{A/Ep)^. The conductance is G ~ A:|, ~ /i, and there- 
fore screening gives a vanishing small contribution to the 
change in the resistance when the splitting is sufficiently 
small, SR^/Ro « {l/48){Ey,/Ep){A/Ep)\ 

Using parameters for Fe, Ni and Co (A™ = 40 nm, 
A^ = 100 nm, A^ = 15 nm, respectively), we find 
R^,/R = 0.0008%, Rn,,/Ro = 0.0001%, R^^./Rq = 0.008% 
respectively. Within the 2-band model, the ballistic do- 
main wall scattering is thus very weak. In first-principles 
band structure calculations these small numbers are en- 
hanced by orders of magnitude due to the (near) degen- 
eracy of the energy bands at the Fermi levelE 



IV. DIFFUSE TRANSPORT 



The DC conductivity of a single-domain ferromagnet 



When the system size is much larger than the mean free 
path, the transport is in the diffuse regime. We assume 
that the electrons are subject to spin-dependent scatter- 
ing, which is modelled by short-range scatters giving rise 
to spin-dependent life-times r+ and t_ for spin-up and 
spin-down states, respectively, which will be treated as 
adjustable parameters. 

We study the current in the ^-direction. The current 
operator transformed by the local spin rotation U{r) in 
(0) is J = W JU = Jo + Jg. The unperturbed current 
operator is 



eh 



m ^ — ^ 



ks'^ks 



ks 



nQ 



(8) 



(9) 



and due to the local gauge transformatio: 

kqss' 

The conductivity is calculated from the Kubo formula. 
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where Uq = N/V is the electron density, N is the number 
of electrons and V is the volume of the system. The 
current-current correlation function n(w) is obtained by 
an analytical continuation (iuji —^uj + iS, 5 —f 0^) from 
the Matsubara correlation function 



n' = --^ 
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where 1/(3 — ksT, ks is the Boltzmann constant, T is 
the temperature and Tt- is the r-ordering operator. We 
will only study the DC conductivity at low temperatures 
by letting uj ^ and T ^ 0. The relevant Feynman 
diagrams to the lowest order in the scattering by the do- 
main wall were identified in Ref. g and are shown in Fig. 
n. Diagram (0) represents the zeroth order Drude con- 
tribution, diagram (1) is due to the correlation of the 
correction to the conductivity operator (0), diagrams (2) 
and (4) are self-energy corrections from the interaction 
Hamiltonian (Eh to the electron Green's function, dia- 
gram (5) is a vertex correction, and diagram (3) is the 
correlation of the change in the current operator (g) and 
the interaction Hamiltonian (H). The electron Green's 
function appearing in the Feynman diagrams in Fig. |l| 
is a configurational average over impurity positions, e.g. 
the retarded Green's function is 



GSsi^) 



huj — e^ + ih/2Ts 



(12) 



The scattering lifetimes of the states at the Fermi en- 
ergy due to the impurity scattering Tji-are assumed to be 
isotropic but may be spin-dependent Jla 



o-o 



vJ: 



V 



dk. 



5(evs - tJ-ojTs 



(n+T+ +7i_r_), 



(13) 
(14) 



where /io is the bulk chemical potential, n+ (n_) is the 
electron density of spin-up (spin-down) states and r+ 
(t_) is the scattering life-time of spin-up (spin-down) 
states. 

There are two contributions to the conductivity which 
to the lowest order in the domain wall scattering are addi- 
tive. Firstly, screening shifts the chemical potential and 
induces a DW conductivity in (^3|). Secondly, the elec- 
trons are directly scattered by the domain wall by the 
interaction term (^ and the gauge transformed current 
operator (^). 

Since the width of domain walls in transition metals is 
much larger than the screening length, electroneutrality 
dictates that the electron density to a good approxima- 
tion is the same in the presence or absence of the DW, 
but the chemical potential differs. This is in contrast to 
the treatment in Refs. ||,0, where the chemical potential 
is assumed to be the same in the presence or absence 
of the DW. The change in the conductivity due to the 
chemical potential shift can be found from (O) setting 
Mo ^ Mo + 5^1-. 
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which to lowest order in the domain wall scattering 
may be added to the DW conductivity. Here Ng — 
mkp / (2it'^Ti?') is the electron density of states at the 
Fermi energy, kp is the spin-dependent Fermi wave vector 
related to the spin-dependent electron density by Us = 
{k^pY / (<oit'^) and we also introduce e|, = h'^ (kj,)'^ / {2m) . 
We proceed by calculating the chemical potential shift 
due to the rotating magnetization. The zeroth order con- 
tribution to the electron density is no — J2ks ^(M^^k)/^' 
where 0{x) is the Heaviside step-function. The Feynman 
diagrams of the contributions to the electron density in 
the lowest order interaction with the domain wall are 
shown in Fig. 0. Diagram (A) is due to the first term in 
(Q) and diagram (B) is due to the second order contribu- 
tion of the second term in (0). Combining the two terms, 
the second order contribution to the electron density is 
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<5(6L-m), (16) 



where 2A(kg) = e^^ — e^_ and k± = k ± (g/2)z. Since 
the DW is much thicker than the Fermi wavelength, we 
disregard the wave- vector dependence {q) on the electron 
states at the Fermi level and introduce the energy param- 
eter for the domain wall, E^ ~ X^o ^^l^9p/(2"^)j ^-9- 



with cos 6' ~ tanh(2/Aiy), E^ — irh^ /{LXwrn), where 
nw = i/L is the "density" of the domain wall. The 
chemical potential shift follows from uq + n2 = n, where 
n is the electron density, 
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The correction in the conductivity (|l^) due to the chem- 
ical potential shift becomes 
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The corrections to the current-current correlation func- 
tion from diagrams (1-5) are 
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where the frequency summations are defined by 
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(19a) 
(19b) 
(19c) 
(19d) 
(19e) 

(20a) 
(20b) 
(20c) 
(20d) 
(20e) 



In the low impurity density limit the energy splitting 
between the bands is larger than the broadening of the 
bands due to the impurity scattering, Ar/Zi ^ 1. The 
frequency sums (|20| ) are evaluated in Appendix ^, where 
the weak wave-vector dependence on the electron states 
at the Fermi level is disregarded consistent with the treat- 
ment of the chemical potential shift above. Carrying out 
the Matsubara frequency sums over the internal energies 
at low temperature, we find that the correction to the 
DC conductivity due to the DW is 
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(21a) 
(21b) 

(21c) 
(21d) 
(21e) 



The contribution to the conductivity from the dia- 
grams (1-5) in Fig. |l] is 
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The DW resistivity can be found from p^ 

where the DW conductivity change due to the rotating 

magnetic field is (7^ — SaQ + J2i=i "'«■ 
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where Sep — J2s sN^ep/ J2s -^s- The first term in (E3|) is 
always positive, but the second and the third terms can 
be negative when the relaxation time of the minority spin 
electrons is longer than the relaxation time for the ma- 
jority spin electrons. However, as will be demonstrated 
below, the domain wall resistivity given by (03) is always 
positive. Our speculation about the possibility of a nega- 
tive domain wall resistance in Ref. IJ is thus not justified 
from Eq. (|2^) only. The result (|23) differs from that ob- 
tained in Ref. ra for spin-independent relaxation times 
Tg = r, where screening was not taken into account, i.e. 
a constant chemical potential and not a constant elec- 
tron density was assumed. The result also differs from 
the calculation in Ref. H based on the Boltzmann equa- 
tion. We believe that this latter discrepancy is because 
in Ref. R screening as well as the effect of the gauge trans- 
formation on the current operator (m are neglected. The 
latter corresponds to the neglect of the diagrams (1) and 
(3) in Fig. nl For sufficiently weak impurity scattering 
and/or a large spin splitting (Ar/h ^ 1) the contribu- 
tion from diagram (1), Eq. (21a), vanishes, so only in 
this limit the omission of this diagram is justified. Fur- 
thermore, the difference in the Fermi wave vectors for 
the spin-up and spin-down electrons was disregarded in 
parts of the calculations in Ref. u\ and an approximation 



was introduced in order to solve the integral equation for 
the Boltzmann equation. Indeed, assuming small spin 
splittings (A/e_F < 1, but Ar/?i > 1) in Eq. (H), we 
obtain 



Po 
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(24) 



which is very similar, but not idential to the result in 
Ref. M. In this limit, the domain wall resistivity increases 
quadratically with the asymmetry in the spin-up and 
spin-down scattering lifetimes as pointed out in Ref. R. 
For larger spin splitting, Eq. ( [2^ ) should be used. 

It is also interesting to study the domain wall resistiv- 
ity when the spin-splitting is large, 2A > fj,, which is the 
case for a half-metallic ferromagnet in which the minor- 
ity spin density of states vanishes iV_ = 0. In this regime 
we find 



Pw 

Po 
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2A 5^2A'^ ^ r_' 



(25) 



The first term in ( |25| ) can be interpreted as additional 
intraband scattering in the majority spin channel due to 
the rotating magnetization and the second term as vir- 
tual transport in the minority spin channel which has 
a negative contribution when t_ > r_|-. The domain 
wall resistivity ( p5| ) is always positive. In the limit of 
large spin-splittings A 3> /^ the domain wall resistivity 
vanishes, since the coupling between the bands becomes 
vanishing small. Note that the present formalism is valid 
only for wide walls, since the domain wall scattering is 
treated as a perturbation. 

Our perturbative result (E3h can be checked against an 
exact calculation for a spin-spiral, dO{z)/dz — oq, a^ ~ 
o-o^q.o with spin-independent life-times Ts -^ r. The de- 
tailed calculation is shown in Appendix |C[ The Hamilo- 
nian is diagonalized in spin space by u± ~ A/±[l,?(l T 
\/l + a^) / a]'^ , where Af± is a normalization constant, 
a — kzUo/p^ and A — Ti p^ /2m. The corresponding 

eigenvalues are E^ — (h^ /2m) f fc^ + a§ =F ^/^fa^~fp^] . 

The (Drude) resistivity can be calculated from the Kubo 



formula. 






(26) 



This is in exact agreement with Eq. (|2^) for Uq = aoSq.o 
when T+ = T_; a good indication of the correctness of 
our perturbation approach. The calculations in Refs. o,M 
disagree with the exact result (26) presumably due to the 
reasons outlined above. 

The result for the domain wall resistivity ( p3| ) can be 
analyzed by introducing kj. = ^J^kp, kp = kp/^/^, 



VVT, 



r/y/rj, where 7 — kp/kp is a mea- 



sure of the polarization of the ferromagnet, 77 = t^/t- 
is a measure of the asymmetry of the scattering life- 
times, kp is the average Fermi wave vector and r is the 



average scattering life-time. The domain wall resistiv- 
ity is proportional to k = E^//AEp. Typically in Fe 
kp ^ 1.7A, Xw ^ 3OOA and n^, r^ 2.5pm~^ giving 
K ~ 10~^, which means that the domain wall scattering 
for symmetric scattering life-times is very weak. How- 
ever, as pointed out in Ref. [7|, the domain wall resistivity 
can become appreciable larger when taking into account 
the life-time asymmetry of the carriers. We show in Fig. 
g the scaled domain wall resistivity Pw/{npo) as a func- 
tion of the asymmetric scattering life-times t^/t- in the 
case of a small spin polarization 7 — 1.01 (solid line), in- 
termediate spin polarization 7 = 1.20 (dashed line) and 
large spin polarization 7 — 10.0 (dotted line). For a 
larger spin-polarization, the domain wall resistivity nat- 
urally becomes asymmetric in the relative difference in 
the scattering lifetimes r_|_ and r_. The domain wall re- 
sisitivity becomes noticable for asymmetric life-times and 
can become of the order Pw/ Po ^ !%• 



V. CONCLUSIONS 

We studied the contribution of domain wall scattering 
on the transport properties of a ferromagnet using an 
effective 2-band model. 

In a diffuse ferromagnet, the domain wall resistivity is 
calculated from the Kubo formula. The domain wall re- 
sistivity is found to be strongly enhanced when the scat- 
tering lifetimes of the majority spins and minority spins 
are different, in agreement with the results in Ref. M. 

In the ballistic regime, we have demonstrated how the 
domain wall scattering creates an effective barrier that 
the electrons must pass. The results from the 2-band 
model give only very small corrections to the resistance 
of the system. 

First-principle band-structure calculations have shown 
that the domain wall resistance can be increased by or- 
ders of magnitude in the ballistic regime.u It would be 
interesting to perform a realistic band-structure calcula- 
tion also for diffuse systems. However, the generalization 
of our 2-band results turns out to be cumbersomeO 
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APPENDIX A: ADIABATIC APPROXIMATION 



In the ballistic regime, it is most convenient to start 



with the Hamiltonian in its first quantized form, whi' 
after the gauge transformation (0) reads H = Hq + V 



tf 



H„ = -T^y^ + Aa,, (Al) 

y = 7^(^ya{z)pz - T—cjya'{z) + -—a^iz) , (A2) 
Zm 2m am 

where a(z) = d0{z)/dz is the gradient of the rotating 
angle of the magnetization and a'{z) — (P9{z)/dz'^. The 
wave function can be written as 



*(r) - 0(p) 



A^iz) 



+ Sk(z) 



(A3) 



where ^(p) — A^^/^ exp(jk||p) is the transverse part of 
the wave function (k = {k^^^kz)), ^k('Z) is the spin-up 
like longitudinal amplitude and Bk(z) is the spin-down 
like longitudinal amplitude. The Schrodinger equation 
then becomes 



and Z — - J dxg'{x), so that the adiabatic solution is 



^k(^) ~ ao(z)exp / dxg'{x) 



(A9) 



Similarly we can find a solution for Bk{z). Disregarding 
tunneling states which only give an exponentially small 
contribution to the conductance, the number of propa- 
gating mo des is determined by the condition Im[g'(a;)] = 
0. From (A8), we see that the number of propagating 
modes is determined by the position where a{z) attains 
its maximum, i.e. the conductance can be calculated as 
for a spin-spiral with a(z) — > Omax- 



APPENDIX B: FREQUENCY SUMMATIONS 

The typical frequency sum to be perfomed is 



7r' = i^X"+'r" + (/^-/), 
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where 



a(z)V4T/ 



(A4) 



(A5) 



ki 



2mEp/Ti — k?,, and Ep is the Fermi energy. The 



off-diagonal terms in (A4) describe the coupling between 
the spin-up and spin-down like states. In the case of a 
spin-spiral, a'{z) = 0, the eigenstates can be found to 
be A\^ — A^eyi\){ikzz), B^ — B^eyip{ikzz), where the 
dispersion of the modes, kz, are determined by 



fci 



kl + aV4 ± Vp^ + q^kl , 



(A6) 



The coupling is weak when the gradient of the spin rota- 
tion gradient is slow compared to the Fermi wave length. 
This permits us to make use pi a multiple scale analysis 
(or adiabatic approximation) Ji3 This analysis is done by 
introducing the small parameter e, so that ijj{z) -^ oj{ez), 
a{z) -^ a{ez) and da/dz — > eda/dz and we introduce 
the newpvariable Z = g{ez)/e, where g{ez) is a scaling 
functionJl3 We expand the longitudinal function Ai^{z) 
and i?k(z) to the lowest order in the small parameter e, 
Ak{z, Z) = al.{Z, z) + 0{e). To the lowest order in the 
small parameter e, the equation to solve is thus 
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We now make the ansatz 0^(2;, Z) — Oj,' (z) exp(iZ) and 
find that the scaling function g is chosen such that 



where X" and Y" are Matsubara Green's functions. 
They can be written in the spectral representation 



X" = 



de Sxje) 
27r icur, — e 



(B2) 



where the spectral function is determined by the retarded 
and the advanced function 



Sxie)^t[X%)~X^ie)] 



(B3) 



Performing the frequency summation in (Bl), we get 

-' = r ^ r pSx{e,)Srie,) x 
/(e2)-/(ei) fiei)-f{e2) 



iuji - (ei - £2) 1(^1 - (£2 - ei) 



(B4) 



The DC-conductivy is obtained by an analytical contin- 
uation 



' lim 



(B5) 



and we consider the limit of zero temperature (T — > 0): 



III 



jRc [^^(0)^^(0) - X"{Q)Y"'{ 
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The product of the two retarded (advanced) Green's 
functions vanishes when integrating over the energy since 
the poles are on the same side of the imaginary plane. 
The sum can then be simplified to 



/ = -Re [X^{0)Y' 



(B8) 



This relation will be used in the following in order to 
calculate the contributions from the diagrams 1-5. 
We use 



Gf,{0)Gf{0) 



G^,{0)-Gm 



(B9) 



where Ss = Ti/(2ts) and obtain in the limit Sg 'ti fJ. the 
contribution from tti to the conductivity 






2i-l ' s I s' 



s{e) + s{n , (BIO) 



where ^^ = e'' — /i is the quasiparticle energy relative to 
the Fermi level. In the case of no spin-splitting (^^ — > ^) 
and {Ss -^ S), the result is /i ~ T^iO- I'^ the limit 
of strong spin-splitting, the result is vanishing small (of 
order Ti/At small). 



/i«0. 
The sum tt^ gives a contribution 

The sum tt^ gives a contribution 



A = 



1 



n {^'Y + 51 {i-^Y + 51 
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In the limit of vanishing spin-splitting and equal life- 
times, the contribution is J3 « T{d/d£)5{£), which agrees 
with the result of 772 as it should. When Ats/Ti ^ 1 
(large spin-splitting) it is 



2r, 



Mn- 
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In the case of no spin-splitting, the sum tt^ gives a 
contribution 



-\ms-' + l5"ios-' 



In the general case, we use 

Gf{0)-G^M 



Gf{0)G^M = 



i {5^s - 5s) - (e- 
For large splitting the result is then 

2r, 
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Finally, the sum tts gives a contribution 
_ h 1 1 

^' ^ 2^WyT5lW^yTW; ■ 

In the case of no spin-splitting the sum is 



\ms-' + l5"i06-' 



For large spin-splitting, we have 
1 



h 



{e'-e-'^y 
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APPENDIX C: SPIN SPIRAL 

The spin-spiral system has a constant gradient of the 
rotating magnetization direction (a^ = 00^9,0)- We per- 
form the local gauge transformation (ph. The trans- 
formed Hamiltonian is 

■^ = 2^ XI ^k ( '^^ ~ '^'P^ + ^okzCTy + -^al j Ck , (CI) 

where c^ is an annihilation operator in the spinor spin- 
space and the exchange splitting A = tYp^/{2m) has 
been introduced. The transformed current operator is 



k 
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The Hamiltonian ( |Cl| ) can be exactly diagonalized, and 
the eigenvalues are 



?± 
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^^=2;^l^' + ro^V^'""+^' 



with the corresponding eigenvectors 

1 



U± =7V: 
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where the parameter a = fc^oo/p^ is introduced and the 
normalization factors are 



JVl = 



2Vl + a2(zpl + Vl + a2) 



(C5) 



The annihilation operators are transformed as c 
(u+, u_) a. In the new basis, the current operator is 



m ^ ^ 



,t ( ^ 2VT+^ 2VT+^ 



"^ k \ 2\/l+Q^ 

The DC conductivity is 



ak (C6) 



2VT+0 



kss' 



where the electron spectral function {A 
the Fermi level is 



Al 



h/T 



[Ei ^lf + {h/ry 



(C7) 



-2IniG'^) at 
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Here we have inserted a phenomenological scattering life- 
time, which is identical for the two eigenstates. Note that 
we cannot treat different scattering life-times for the mi- 
nority and majority states in the bulk ferromagnet with 
the method outlined in this appendix, since the life-times 
appearing in (C7) are the life-times for the exact eigen- 
states in the spin-spiral. In order to determine the rela- 
tion between the different lifetimes, the general met hod 
described above in our paper should be used. From (C7) 
and ( |C8| ) it can be seen that he off-diagonal terms in the 
conductivity, A+A- are in the order l/ir/S./fif' smaller 
than the diagonal terms. We further assume that the 
scattering by the domain wall is weak, i.e. a§ ^ p^ and 
cn)kp <C p^ and expand the result for the conductivity to 
the second order in ag. The conductivity becomes 



^2 _ ^ 
67r^m -^-^ ^ ^ 
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where fi k /{2m) — ii and E^ = h a^/{2m) are used. 
The conductivity should be related to the electron den- 
sity by eliminating any reference to the chemical poten- 
tial which may change in the presence of the domain wall, 
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Inserting ( |C10|) into (|C9| ), the conductivity can therefore 
be written as 
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cr = cro 1 
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where aQ = e^(n_|_ -I- n_)T /m is the Drude conductivity. 
The domain wall resistiviy, p^ — —ba^ja^ is thus 






where po — I/cq. 
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FIG. 1. Feynman diagrams of the contributions to the con- 
ductivity to the lowest order in the domain wall scattering. 
Solid lines indicate the electron Green's function and the 
dashed lines the interaction with the domain wall. The ver- 
tex X arises from the unperturbed current operator (Bl), the 
vertex o is due to the gauge transformation on the current 
operator (19), the vertex o is due to the first term in the inter- 
action Hamiltonian (0) and the vertex □ is due to the second 
term in the interaction Hamiltonian (m. 



t also at Philips Research Laboratories, Prof. Holstlaan 4, 
5656 AA Eindhoven, The Netherlands. 



FIG. 2. Feynman diagrams of the contributions to the elec- 
tron density to the lowest order in the domain wall scatter- 
ing. Solid lines indicate the electron Green's function and the 
dashed lines the interaction with the domain wall. The vertex 
o is due to the first term in the interaction Hamiltonian (M) 
and the vertex □ is due to the second term in the interaction 
Hamiltonian (0). 



FIG. 3. The relative change in the resistiviy, p^/{iipo), due 
to the scattering by the domain waU as a function of the asym- 
metry in the scattering hfe-times t+/t-. The sohd line is for 
7 = kp/kp = 1.01, the dashed line is for 7 — kp/kp = 1.20 
and the dotted line is for 7 — kp/kp — 10.0. 
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